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It is studied the scattering of magnons by the 2d topological Belavin-Polyakov soliton in isotropic 
ferromagnet. Analytical solutions of the scattering problem are constructed: (i) exactly for any 
magnon wave vectors for the partial wave with the azimuthal number m — 1 (translational mode), 
and (ii) in the long- and short-wave limits for the rest modes. The magnon mode frequencies are 
found for the finite size magnets. An effective equation of the soliton motion is constructed. The 
magnon density of states, connected with the soliton-magnon interaction, is found in a long-wave 
approximation . 
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Nonlinear topologically nontrivial excitations (soli- 
tons) are well-known to play a special role in a low- 
dimensional magnetic systems. For example, the pres- 
ence of vortices in 2d easy-plane (EP) magnets gives rise 
to Berezinskii-Kosterlitz-Thouless phase transition 0. 
Kinks in Id magnets and localized Belavin-Polyakov soli- 
tons (BP-solitons m ) in 2d isotropic magnets are respon- 
sible for the destruction of the long-range order at finite 
temperature. This can be explained within the scope 
of so-called soliton phenomenology, where the magnet 
can be described as a two-component gas of elementary 
excitations: solitons and magnons. Such approach was 
developed for Id magnets see also Rcfs. |^J^. The 
soliton signature in dynamical response functions can be 
observed experimentally. Translational motion of soli- 
tons leads to the so-called soliton central peak. An- 
other possibility to detect the soliton signature is to look 
for magnon modes, localized on the soliton (local modes, 
LM), observed in For Id magnets such scattering 
causes the change of the magnon density of states which 
is necessary for self-consistent calculation of temperature 
dependence of the soliton density 1^ . 

For 2d the concept of soliton-magnon gas has been ex- 
tended, e.g. to describe EP magnets 0, and to explain 
the EPR line-width in easy-axial magnets How- 
ever the general behaviour of the 2d soliton dynamics is 
not clear at present. In particular, the form of inertial 
terms in the dynamical equations for the soliton centre is 
unknown; the soliton density have not been calculated, 
but been used as input parameter. 

The problem of the soliton dynamics or the problem of 
LM existing are intimately connected with the soliton- 
magnon scattering. For example, using numerical data 
for the scattering amplitude the non-Newtonian effec- 
tive equation of motion of the magnetic vortex was con- 
structed ||l^. For the 2d EP antiferromagnet (AFM), 
finite-frequency truly localized internal mode was pre- 
dicted ||ll|. For all mentioned papers, an analysis of 2d 
solitons was carried out numerically. It becomes espe- 



cially important to analyze models for which analytical 
results can be obtained. We are aware of only the exact 
solution for BP-soliton ^ , which describes a topological 
soliton in isotropic magnets. 

In this Letter we have investigated the soliton-magnon 
scattering for the BP-soliton in the isotropic ferromag- 
net (FM). An exact analytical solution of the scattering 
problem is constructed. An effective equation of the soli- 
ton motion is found for the finite size FM. We also ana- 
lyzed the long-wave asymptotics of the magnon density 
of states in the presence of BP-soliton. 

Solitons and Magnons. The dynamics of the FM 
follows the Landau-Lifshitz equations for the normal- 
ized magnetization m JT2| . In angular variables, nix + 
imy = siii9exp{i(j)), these equations correspond to the 
Lagrangian L = {A/2) f d^xC, 
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where the A = JS^, J is the exchange integral, and S 
the atomic spin. The magnon solutions have a form 9 — 
const <^ 1, (j) — k ■ r — uj{k)t with the dispersion law 
Lo{k) — Dk^, so the constant D is the spin-wave stiffness. 

The simplest static nonlinear excitations in 2d case are 
the BP-solitons B: 
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Here r, x £^re the polar coordinates in the magnet plane, 
the integer v determines 7r2 topological charge of the soli- 
ton, and R , (/Sq a-re arbitrary parameters. The energy of 
the soliton (||) Eq — A'kA\v\ does not depend on the soli- 
ton's radius Rq due to the conformal invariance of the 
static model Apparently, such invariance breaks for 
FM in the dynamical case. 

The scattering problem. The equations for linear os- 
cillation on the soliton background, 9 = 9o{r) + t?(r, x, t) 
and (j) — + + (sm 9o)~^iJ,{r, x, t) can be represented 
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in the form of the single equation for the complex quan- 
tity ip = "d + ifi, 



_V2 + ^cos20o)V'-5-cos0o^ = TT^. (3) 



2iiy 



dip i dip 
'9^ " D~dt 



The solution of Eq. (j^) has the form of a superposition of 
cylindrical waves ip — X]m=-oo /me*™-^"*'^*, m is the az- 
imuthal quantum number. The partial waves fm satisfy 
the 2d-radial Schrodinger eigenvalue problem (EVP), 

Hf„, - [kRffrn, H^~Vl + UM (4a) 
U„i{x) ^ (l/x^) {m^ + 2mi^cos9o + cos 26*0) . (4b) 

Far from the soliton, r ^ R, the function Um ~ /x^ — 
4z/(to + 2i/)(1/x)2(''+i), so 



/„ cx J|p| (z) + (z) + Av{m + 2v) {kRY'' F\p\ (z), 
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where z = kr, J|p| and are Bessel and Neumann 
functions, respectively, p — v + m. The quantity 
cr^ = — tan 6'^ determines the amplitude of the soliton- 
magnon scattering, the S'-matrix can be written as — 
exp(2i(5J^). The correction F, which is small at kR ^ 1, 
has a power decay F cx z~^^^ at kr ^ 1, but its contribu- 
tion is important to calculate the scattering amplitude. 
The explicit form of F can be found after long but simple 
exercises with the unhomogeneous Bessel equation. 

In the case of w = there is an exact solution of 
EVP (P due to the restoration of the conformal invari- 
ance of the model : 
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(we shall discuss the case > 0; for analyzing of < 0, it 
is sufficient to replace m by — m). It should be noted that 
solutions (^ are finite and have no singularities at r ^ 0, 
as well as at r — > 00 when — < m < This at once 
leads to the existence of 2|i^| LM with zero frequencies 
(zero LMs). The physical meaning of two such modes 
is obvious: the translational mode f^n=i describes the 
displacements of a soliton as a whole; the mode /^Iq 
corresponds to the change of (fio and soliton radius R. 
Note that LMs are limit points for magnon modes of the 
continuous spectrum when k (unlike Id magnet, in 
which LMs are separated from the continuous spectrum 
by the gap [|). 

For m > V, the solution fm^ has a singularity in 
the origin, but in this case the second independent so- 
lution of EVP (|) with fc = , regular at r ^ 0, 
/W = a;™sin6'o [a;2'^/(m + i^) + 2/m + a;-2''/(m-zy)] , 
can be used. Thus, for = 0, we can construct at least 
one solution fm^ or fin of the EVP (^ that does not 
have a singularity in the origin. 



The existence of the exact solution fm^ permits to 
present the Schrodinger operator H in the factorized form 
H = iti^ i/^o) ^ 0, where 



A = 
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Introducing the function 
rewritten as 

ngm^{kRfg,n, n = AA^ 



Afm, the EVP m can be 
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The function /,„ can be found via g„i as /,„ = 
{kR)~^A'gm . Using such transformation one can find 
an exact solution for the translational mode. Indeed, in 
this case Ui{x) = i^^/x^, so gi = —J^{kr) and 



h = J. 
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The existence of the exact solution for any k is the unique 
property of the model ([|). It demonstrates the fact, men- 
tioned above (||), that the deviation /i from the asymp- 
tote Jjy+i is specified by power decaying apart from the 
soliton (while the deviations far from the vortex core are 
decaying exponentially [0,|o|). Again, the EVP (||) is 
more convenient for numerical analysis because the "po- 
tential" Um is repulsive. 

In order to describe the magnon scattering in the long- 
wave approximation we use the fact that exact zero func- 
tions /m-* or /m-* describe correctly the solution in the 
region r <^ 1/k. For small but finite values of k wc make 
the ansatz 
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where is one of the exact zero solutions, seen above. 
On the other hand we can use the scattering approxi- 
mation with account of the correction F (^). Therefore, 
over a wide range of values of r, for R ^ r <^ we 
can use the ansatz (|l0|): on one hand, starting from the 
exact zero solution, and on the other hand explaining the 
results in term of the scattering problem. Awkward cal- 
culations lead to long-wave asymptotics of the scattering 
amplitude cr!^^^^ for the soliton with u — 1 having mini- 
mal energy. It was found that the scattering intensity is 
maximal for the LM with m = 0, 
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, fciZ<l. 



(11a) 



The same dependence occurs for solitons with 7^ 1 at 
TO = -~i> + 1. 

The second LM (translational mode, m = 1) does not 
scatter in accordance with Eq. (M). For the quasi-LM 
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(to = —1) we can write (T_i(fc) = 7r(fci?)^ In (1/fci?). In 
the range of absence of LM, we can restore the general 
dependence 



O-m(fc) 
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fci?< 1. 
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In the short-wave hmit k ^ \m\/R we can use WKB 
approximation for all r far from the turnover point, 
|TO|/fc <^ r, and at the same time the asymptotic gm oc 
Jli'-m+ilicr) at r ^ i?. Thus at a wide range of values 
of r, \m\/k <C r ^ i?, both solutions are valid. Their 
comparison gives: 
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One can see that cr„i tends to zero for both long- and 
short-wave limits, but with opposite signs. Due to this, 
the values of (5J^ at fc ^ cxd and k ^ Q differ by ±7r, the 
sings "+" and "— " correspond to to > 1 and to < 0, 
respectively. Therefore, the scattering amplitude has a 
pole at some finite k = kp, the value of kp is growing 
as |to| with the growth of |to|, which is confirmed by 
the numerical calculations. Note that for the magnon 
scattering by Id soliton in a number of models the values 
of (5 at fc — ^ CX3 and A; — > differ by — tt, that causes the 
decrease of the total number of magnon states by one . 

Finite size magnet. The investigation of the magnon 
eigenmodes on the soliton in the finite size magnet can be 
used for the analytical description of the direct numerical 
simulation data of the soliton motion. In this way the 
non-Newtonian vortex dynamics was explained in ||ic| ] 
for EP FMs. On the other hand such calculation can 
be used to describe eigenmodes for small FM particles in 
so-called vortex state |0. 

We consider magnon modes in the circular system 
of radius L with the soliton in the centre. Differ- 
ent boundary conditions (BC), both fixed (Dirichlet 
BC, DBG, fm\^^^ = 0) and free (Neumann BC, NBC 
{dfm/dr)\^_j^ = 0) will be discussed. Without soliton 
the magnon spectrum in such a magnet is discrete, e. g. 
for DBC the eigenvalues km,i = jm,i/L, where is 
the i— th zero of J,„. If the soliton is present, this is 
true again for some modes with small scattering <Tm ^ 1 
km,i = Zm,i/L, where z^a is the i—th root of the 
equation J|p| (z) + cri^p| (z) = 0. 

However, in the LM case < m < v, the system has 
higher symmetry due to the restoration of the confor- 
mal invariance. Naturally, in the finite size magnet case 
there are quasi-Goldstone modes (qCMs) with anoma- 
lously small frequencies, i. e. kL <^ 1. In particular, the 
qGM appears for the translational motion (|to| = 1) of 
the vortex in the EP FM. The existence of such a mode 
is determined by the scattering only [|ll],0 . 

We discuss the most important case of qCM with to = 
1, which describes the translational motion of the soliton 



and could be used for the construction of the dynamical 
equations of the soliton ||l^,|l^ . Using the exact solution 
(^ in the region kR <C 1 one can find the frequency of 
the translational qGM in the form 



ujQ = ±8DR^L-'^ 



(12) 



signs "— " and correspond to DBG and NBG, re- 
spectively. The next translational mode has a frequency 
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or wi 



Dj'^^L ^, respectively. The 



highest modes are separated from these ones by the finite 
gap, so the picture of doublets, which is characteristic for 
the vortex-magnon scattering, is absent [|l0|. Thus it is 
a belief that the account of these two modes gives the 
adequate description of soliton motion. 

To describe the soliton dynamics with two character- 
istic frequencies we can use the 2'"* order equation 
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where M and G are effective mass and the constant of 
gyroforce, G = AttvA/D is determined by the topology 
only [|l6|,|j , F is an image-force acting on soliton due to 
the boundary. Such a form of F can be explained by 
taken into account that the magnetic vortices interact 
as 2d charges, and the BP-soliton with v ^ 1 can be 
interpreted as a vortex dipole. 

Assuming ujq <C wi, it is easy to compare wq, i^i with 
the frequencies of the Eq. (|l^ and extract the unknown 
value of the effective mass (this method was used for the 
vortices in EP magnets with the using of the numerical 
values of a;o,^^i as an input [psf). Indeed, cjq ~ — a/G 
coincides with the frequency of qGM. Comparing uji with 
the value ui ~ —G/M we obtain 



Mn = 



An A f L 



D2 



Mn = 



An A ( L 



D2 



(14) 



for DBG and NBG, respectively. Therefore M is unlocal 
as the coefficient G3 in 3'''^-order equations of motion for 
the FM 0. Note that the dependence M oc agree 
with the result M o^XjK^K \s the anisotropy constant, 
for the easy-axis FM |]l^, if the characteristical length 
Ao = \/ A/K is changed by L. 

The analysis shows that qGMs appear for all higher v 
in the region of the existence of LMs, —v < m < v. Their 
frequencies are proprotional to (Z?/L^)(i?/L)^l''+'"^^l or 
[D/L^ ln(L/i?)]. 

Magnon density of states. The main point of the 
Id soliton (kink) phenomenology |^ is a change in the 
magnon density of states due to the presence of the kink, 
Pid{k) — {l/2n) dS{k)/dk, where S{k) is a phase shift 
by kink-magnon interaction. It causes the decrease of 
the total number of magnon states, so the magnons free 
energy changes accordingly. 



3 



Let us transfer this approach to the 2d case. Using the 
magnon density of states the sohton phenomenology of 
2d magnets could be constructed. In particular, the soli- 
ton density could be calculated. The free 2d magnon has 
an expansion in terms of cylindrical waves, Jm{kr)e^™-^ , 
with the quantized angular part. Hence only the radial 
part J{kr) must be quantized. In the circular geometry 
with the radius L the eigenvalues fc„ satisfy the condi- 
tions knL — jm,n for the DBC. For n ^ 1 the standard 
condition fc„L « irn takes place. However, for \m\ ^ 1 
the 1"* zero jm.i ~ so the set of allowed values of m 
must be limited by the condition |m| < kL. Accordingly, 
the 2d density of states is 



Pik) 



kL 



m— — kL 



dk ■ 
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For small k, the function p{k) diverges due to the 
mode with m — Q. This contribution is dominant for 
low temperatures. Using ( |lla| ) one can obtain p{k) « 
~ {2k) ^\n~^(kR). However, our numerical calculations 
done for largest k shows that the series (^) is alternating 
and the total number of magnon states does not decrease 
as for Id. There is only a redistribution of the magnon 
modes between the states with the opposite m. 

Conclusion. We have constructed an exact analytical 
solution of the soliton-magnon scattering problem for the 
isotropic FM. We used these results for different proper- 
ties of solitons and LMs. In particular, LM frequencies 
have been calculated for the finite size FM. Eigenmodes 
with anomalous small frequencies appear in the small FM 
particle having cylindrical geometry with the soliton in 
the centre. The effective equation of the soliton motion 
have been constructed. It was analyzed the magnon den- 
sity of states in the presence of the soliton. 

All the results for the magnon modes (obtained for the 
FM) can be expanded to the Lorentz-invariant nonlin- 
ear tr-model. Lagrangians for this model results from 
(1) after replacing of the gyroscopical term with \/ D by 
{l/c^){dl/dt)^, where c is the characteristic speed. For 
magnets, this model can be used to describe the AFM dy- 
namics, where the dispersion law LUpMik) — Dk^ must 
be replaced by loafm — ick. Two lowest modes with 
m — 1 having the frequencies luq — ±2^/2ci?i-2 de- 
scribe two independent modes of translational dynam- 
ics of AFM BP-soliton with G = 0, the effective mass 
Mafm — Eq/c? under the action of image soliton force. 
The magnon density of states p{k) has the same form as 
for the FM. 

We may use these results for the Euclidean version 
of the cr-model, important for the quantum chains with 
the AFM interaction. Properties of such chains are con- 
nected with the instantons of the cr-model. Both instan- 
tons with the BP-soliton structure [ p^ and merons with 
half-integer topological charge ||l^] are discussed. The 



knowledge of the total set of eigenvalues on the instan- 
ton background, especially zero modes, is necessary to 
calculate the fluctuation determinant. 
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